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ABSTRACT 

We present observations of Ti II A3384 absorption towards 15 distant stars in the Galactic 
halo and the Magellanic Clouds. These new data extend existing surveys of the distribution of 
Ti + to larger distances from the plane of the Galaxy than sampled previously, allowing the scale 
height of titanium to be determined for the first time. We find h Ti + = 1.5 ± 0.2 kpc, a value 
which although greater than those of other tracers of neutral gas, is not as large as had been 
suspected. We interpret the extended distribution of Ti + as an indication that its severe depletion 
in interstellar clouds in the disk is reduced at the lower densities prevailing in the halo. The data 
are consistent with a simple power-law dependence of the Ti abundance on the ambient density, 
with exponent k ~ —1. If the model is correct, it implies that refractory elements like Ti are fully 
returned to the gas phase at distances beyond ~ 1 kpc from the plane of the Galaxy. 
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1. INTRODUCTION 

The determination of the extent, distribution, and kinematics of the interstellar medium away 
from the Galactic plane is important both for understanding the physical processes which drive 
the gas to large distances from the disk (e.g. Houck & Bregman 1990) and for interpreting the 
absorption produced by more distant galactic haloes in the spectra of background QSOs (e.g. 
Bowen, Blades, Sz Pettini 1995). This topic has been an active area of research over the last fifteen 
years, but the overall picture which has emerged is still incomplete. In particular, it is not yet 
clear how to bring together the information provided by different tracers of the ISM, which can 
exhibit significant differences in their distributions. 

Highly ionized gas, discovered from C IV and Si IV absorption with IUE (Savage & de Boer 
1979, Pettini & West 1982), has been shown by HST observations of extragalactic sources to have 
an exponential scale height h ~ 5 kpc (Savage et al. 1993; see also the review by Sembach & 
Savage 1992). The patchiness of this gas, however, is such that the observations do not provide 
tight constraints on its large-scale distribution in the Galaxy; a spherical halo fits the data just as 
adequately as the more generally assumed exponential fall-off with distance from the plane, | z \. 

Even the distribution of neutral hydrogen is open to different interpretations, despite the 
wealth of data collected with the IUE satellite and recently re-analyzed by Diplas Sz Savage 
(1994a). By considering the values of hydrogen column density, iV(H°), for 374 sight-lines through 
the Galaxy, these authors deduced a scale height h H o ~ 0.2 kpc, although the observations are 
equally well represented by the sum of two exponential distributions, with hi = 0.12 kpc and 
hi = 1.2 kpc respectively. That significant amounts of H° do exist at | z | distances greater than 
1 kpc is demonstrated by the comparison between 7V(H°) measures from Lyman a absorption 
towards halo stars and from 21 cm emission, the latter sampling the full extent of the gas in a 
given direction, including material beyond the stars (Lockman, Hobbs, Sz Shull 1986; Danly et al. 
1992). 

In disk clouds, refractory elements are mostly removed from the gas phase and incorporated 
into dust grains. Edgar Sz Savage (1989) showed that the scale heights of species such as Ti + , 
Ca + and Fe + are significantly larger than that of H°, in broad agreement with earlier findings 
that the depletions of these elements exhibit a dependence on < n(H°)>, the average gas volume 
density along a line of sight (e.g. Savage Sz Bohlin 1979; Phillips, Pettini, Sz Gondhalekar 
1984). Presumably, as n(H°) decreases with distance from the plane of the Galaxy, the reduced 
depletions of Ti, Ca and Fe compensate for the overall decrease in gas density, leading to shallower 
profiles in | ^ | for these species than for neutral hydrogen. Spitzer (1985) and Jenkins, Savage, Sz 
Spitzer (1986) proposed that the correlation of depletion with mean line-of-sight density can be 
understood in terms of a simple model of the ISM consisting of diffuse cold clouds and of a warm 
intercloud gas. In this picture the refractory elements are more severely depleted in the clouds 
than in the intercloud medium; the observed correlations follow from the fact that < n(H°)> is a 
rough measure of the relative proportions of the two phases along a given line of sight. The two 
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components may well have different scale heights, with the lower density intercloud gas conceivably 
extending further into the halo than the cold clouds. Species such as Ca + and Ti + are then useful 
tracers of a diffuse interstellar medium which may be found out to large distances from the disk. 

Ti + is especially favoured in this context. Its ionization potential (13.58 eV) is so close to 
that of H° that concerns about unobserved ion stages are minimized and the Ti abundance in the 
neutral gas can be deduced directly. Furthermore, the strongest transition of Ti II, at A3383.768 
with an /-value of 0.3401 (Morton 1991), is normally unsaturated, allowing a reliable measurement 
of the column density from the line equivalent width. These properties contrast with those of other 
interstellar species, such as Ca + and Na°, accessible to ground-based telescopes. On the other 
hand the line occurs in the near-ultraviolet, where atmospheric extinction is high and in general 
the sensitivity of optical instrumentation is low. Therefore, observations of Ti II absorption in the 
Galactic halo have until recently been limited to the brightest stars and relatively little data are 
available for distances beyond 1 kpc. 

In particular, there are hardly any published measurements of Ti II in the spectra of 
extragalactic sources, allowing the total column density perpendicular to the Galactic plane, 
7Vo(Ti + ), to be deduced. Consequently, previous attempts to determine the scale height of 
Ti + have been inconclusive. Edgar & Savage (1989) reported a lower limit h Ti + > 2 kpc, 
while h Ca + = 1 kpc. By adding new observations of 25 sight-lines to earlier compilations in a 
comprehensive study of Ti II and Ca II in the halo, Albert et al. (1993) confirmed the scale height 
of Ca + , but were still unable to see a turn-over in the vertical distribution of Ti + . 

In this paper we remedy this situation by reporting observations of Ti II absorption in the 
spectra of 15 stars located beyond 1 kpc from the plane of the Galaxy and including four stars in 
the Magellanic Clouds. These data extend previous studies of the distribution of Ti + in the halo 
to a largely unexplored region and, as we discuss below, show that the scale height of Ti + is in 
fact h Ti + ~ 1.5 kpc, greater than that of Ca + , but not by a large factor. 
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2. OBSERVATIONS AND DATA REDUCTION 

Table 1 lists relevant properties of the 15 stars observed in this study, including their spectral 
types and distances from the Sun. Nominally, the distance estimates are accurate to about ±25%, 
but the question of whether these stars have been correctly classified is still open. Although 
the presence of young B stars several kpc from the plane is surprising, the work of Keenan and 
collaborators in particular has shown that many of these stars have temperatures, gravities and 
chemical abundances indistinguishable from their disk counterparts and are therefore probably 
bona fide Population I stars (Tobin 1991 and Keenan 1992 give recent discussions of this problem). 
In the present work we have adopted, whenever possible, distances deduced by fine analyses of 
high dispersion stellar spectra, such as those published by the Belfast group (see Little et al. 1994 
for a description of the method). For the stars in Magellanic Clouds we have adopted distances of 
50 kpc to the LMC and 57 kpc to the SMC from the review by Feast (1991). 

The data reported here were secured over a period of 5 years, from 1989 to 1994. All but one 
of the stars were observed with the University College London echelle spectrograph (UCLES) at 
the coude focus of the 3.9 m Anglo- Australian Telescope at Siding Spring Observatory, Australia. 
The observations were mostly carried out in periods of poor observing conditions, as back-up to 
a more demanding programme of faint-object spectroscopy which has been described elsewhere 
(Pettini et al. 1990). The spectrum of one star, HDE 233622, was obtained with the Utrecht 
echelle spectrograph (UES) at the Nasmyth focus of the 4.2 m William Herschel Telescope on 
La Palma, Canary Islands; this instrument is very similar to UCLES. With both UCLES and 
UES we used the 79 grooves mm -1 echelle gratings, which result in an inter-order separation in 
excess of 20 arcsec at the wavelength of interest here, thereby providing adequate sky signal for 
accurate background subtraction. All the spectra — except those of HD 5980 and HD 269333 — 
were recorded with the IPCS, a photon counting detector which, because of its low noise, high 
quantum efficiency in the ultraviolet and small pixels, is well suited to the present programme. 
HD 5980 and HD 269333 were observed in 1994 with a 1024 x 1024 Tektronix CCD. The slit width 
was set to 0.8 arcsec for the IPCS observations and increased to 1 arcsec when using the CCD 
which has larger pixels. 

The two-dimensional echelle spectra were reduced with standard procedures using the 
STARLINK software packages FIGARO, DIPSO and APIG; a detailed description of the data 
reduction steps is given by Pettini et al. (1990). In general, several exposures were obtained 
for each star in Table 1. Individual spectra, wavelength calibrated and sky-subtracted, were 
rebinned to a common Local Standard of Rest (LSR) wavelength scale before being added 
together. The profiles of emission lines from the Th-Ar hollow cathode lamp used for wavelength 
comparison show the resolution of the final spectra to be FWHM = 0.060 ± 0.003 A, sampled with 
~ 2.4 IPCS pixels; at the wavelength of Ti II A3383.768, the corresponding velocity resolution is 
5.3 ± 0.3 km s _1 . The resolution of the two CCD spectra is 25% lower (FWHM = 6.6 km s" 1 ). 



TABLE 1 

The Stars Observed 



Star 


/ b 


R.A. 


Dec. 


Spectral Type 


V 


B-V E(B-V) d z vsmi Refs. 


Date 






(1950) 


(1950) 






(kpc) (kpc) (kms -1 ) 


Observed 



HD 5980 


302.1 


-44.9 


00 


57 


22.9 


-72 26 36 


WP + OB + Neb: 


11.81 


-0 


.22 





,09 


57" 


-40" 




1,2,3 


21 Aug 1994 


HD 121968 


334.0 


55.8 


13 


56 


15.8 


-02 40 19.9 


Bl V 


10.31 


-0 


.19 





,07 


5.2 


4.3 


>150 


4 


12 k 13 Aug 1991 


HD 125924 


338.2 


48.3 


14 


20 


03.77 


-08 01 15.9 


B2 IV 


9.68 


-0 


.19 





,05 


3.2 


2.4 


105 


4 


11 Aug 1991 


HD 158243 


337.6 


-10.6 


17 


27 


03.97 


-53 26 28.1 


Bl lab 


8.15 





.00 





,19 


6.5 


-1.2 




5 


14 Aug 1991 


HD 163522 


349.6 


-9.1 


17 


55 


00.23 


-42 28 55.9 


Bl la 


8.44 





.00 





,19 


9.3 


-1.5 


~ 50 


5,6 


3 Jun 1990 


HD 179407 


24.2 


-10.4 


19 


10 


04.94 


-12 40 04.9 


B0.5 lb 


9.30 





.06 





,28 


5.4 


-1.0 


156 


7 


10 k 13 Aug 1991 


HD 195455 


20.3 


-32.1 


20 


29 


17.22 


-24 24 15.7 


B0.5 III 


9.20 


-0 


.18 





.08 


3.1 


-1.6 


200 


7 


14 Aug 1991 


HD 206144 


34.8 


-45.1 


21 


37 


47.31 


-17 49 39.2 


Bl III 


9.41 


-0 


.15 





,11 


5.0 


-3.6 


64 


7 


10 Aug 1991 


HD 209684 


43.5 


-49.1 


22 


02 


49.36 


-14 00 49.8 


B3 V 


9.86 


-0 


.20 





,00 


1.7 


-1.3 




4 


26 Aug 1990 k 14 Aug 1991 


HD 220787 


67.8 


-64.4 


23 


24 


09.74 


-11 18 27.2 


B3 III 


8.30 


-0 


.17 





,03 


1.4 


-1.3 


27 


7 


21 Jul 1989 


HDE 233622 


168.2 


44.2 


09 


18 


08.26 


50 18 43.9 


B1.5 


9.97 


-0 


.23 





,02 


3.5 


2.4 


240 


4 


5 Dec 1992 


HD 269006 


282.8 


-34.2 


05 


02 


44.2 


-71 24 15 


Peculiar 


9.83 










50* 


-26* 




1,3 


14 Aug 1991 


HD 269333 


279.9 


-33.4 


05 


18 


39.2 


-69 14 45 


W + B1:I 


11.21 


-0 


,09 





,11 


50* 


-26* 




1,8,3 


21 Aug 1994 


JL 212 


303.6 


-61.0 


00 


46 


48.35 


-56 22 09.4 


B2 


10.32 


-0 


,20 





,04 


3.6 


-3.2 


219 


7 


26 Aug 1990 


R 136 


279.4 


-31.7 


05 


39 


04.4 


-69 07 35 


OB(n)+ 
WN5-A(B) 


10.77 





.14 





.34 


50* 


-26* 




9,3 


24 Feb 1990 



" In the Small Magellanic Cloud 
4 In the Large Magellanic Cloud 

REFERENCES — (1) Songaila, Blades, Hu k Cowie 1986; (2) Fitzpatrick 1985; (3) Feast 1991; (4) Danly 1989; (5) Savage k Massa 1985; 
(6) Savage, Sembach k Massa 1990; (7) Keenan, Dufton k McKeith 1982; (8) Walborn 1977; (9) Pettini k D'Odorico 1986. 
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TABLE 2 

Equivalent Widths and Column Densities 



otar 


L 




z 


Q /TVT 
b/lN 


TT/ /I \ 

w xyooo^t) 


AT ( r V\~ 


) 


7V(H ) 


l 






(kpc) 




(mA) 


i 1 nil 
(W cm 


J 


f 1 n20 

cm 


) 


nU 5980 


A A a 

— 44. y 


A n 
—40 


1 o 
12 


A 1 l 1 1 

41 ± 11 


1 o n i 

13.0 ± 


o c 
3.5 


1 n 
1.9 




un i oi nco 


K K O 

55.0 


A O 

4.3 


O /l 

34 


OO 1 o 

3o ± 2 


1 o n _i_ 
12. U ± 


n o 
0.8 


A (\ 1 

4.U ± 


I). 8 


HD 125924 


/lOO 

4o.3 


O A 

ZA 


OO 

2o 


/I C _1_ K 

4o ± 5 


1 A o I 

14.o ± 


1.0 


/i i 
4.b ± 


u.y 


Tjr\ -i roO/1 O 

HD 15o24o 


— 10.6 


1 o 
— 1.2 


on 
30 


m i /i 

yi ± 4 


on 7 _i_ 
30. < ± 


1 /I 

1.4 


-ion i 

12. y ± 


o.l) 


un 1 coo 


n 1 
—9.1 


— 1.5 


57 


178 ± 3 


b0.4 ± 


O 7 

2.7 


-1 o O 1 

12.3 ± 


Z.J 


HD 179407 


1 n /i 
— 10.4 


1 n 
— 1.0 


O E 

35 


OO 1 A 

83 ± 4 


27.2 ± 


1 o 

1.3 


-ion i 

12.9 ± 


o o 
J.z 


MD 195455 


OO 1 

—32.1 


— 1.0 


O E 

35 


49 ± 3 


1 F O I 

15.3 ± 


n n 
0.9 


A T 1 

4.7 ± 


1 1 

1.1 


un oAt 1 1/1/1 
HD 20bl44 


/i k 1 
—45.1 


o c 
—3.0 


31 


c 1 i /i 
51 ± 4 


lo.o ± 


1 o 

1.3 


O o 1 

3.2 ± 


1 n 
l.U 


HD 209684 


-49.1 


-1.3 


33 


51 ± 4 


16.2 ± 


1.3 






HD 220787 


-64.4 


-1.3 


40 


26 ± 3 


8.0 ± 


0.9 






HDE 233622 


44.2 


2.4 


55 


27 ± 3 


8.1 ± 


0.9 






HD 269006 


-34.2 


-26 


12 


55 ± 12 


17.3 ± 


3.8 






HD 269333 


-33.4 


-26 


20 


28 ± 4 


9.3 ± 


1.4 


4.3 




JL 212 


-61.0 


-3.2 


22 


42 ± 5 


13.0 ± 


1.5 






R 136 


-31.7 


-26 


26 


62 ± 7 


21.5 ± 


2.4 


3.2 





HI data from Lyman a measurements by Diplas & Savage (1993a) except for HD 5980, HD 269333 and 
R 136 which are from the 21cm observations by McGee, Newton & Morton 1983. 

Portions of the echelle spectra encompassing the interstellar Ti II line are reproduced on 
an expanded vertical scale in Figure 1; the spectra have been normalized to the local stellar 
continuum. As can be seen from the Figure, although interstellar Ti II absorption is detected in 
every case, the signal-to-noise ratio varies significantly from star to star. Values of S/N, deduced 
from the rms deviation of the data points from the fitted stellar continuum, are listed in column 
4 of Table 2. Typically, S/N ps 30; the corresponding 3cr detection limit for the equivalent width 
of an unresolved absorption line is W\(Sa) ~ 6 mA. The measured equivalent widths of the Ti II 
lines are listed in column 5 of Table 2, together with la errors which reflect only the counting 
statistics (and do not include uncertainties in the continuum placement). From Figure 1 it can 
be seen that in the case of R 136, and possibly HD 269006, we detect Ti II absorption at LMC 
velocities of between vlsr — 200 and 300 km s _1 ; in HD 5980 there is a 3.5cr absorption feature 
at SMC velocities (ps 140 km s" 1 ). 



HD 269333 



JL 212 



1.2 



1.0 



0.8 




1.2 



1.0 



0.8 




-400 -200 



200 400 



-400 -200 



200 400 



1.2 



1.0 



0.8 



R 136 




-400 



400 



Fig. 1. — Normalized profiles of Ti II A3383.768 absorption lines towards 15 distant stars, including one in 
the Small Magellanic Cloud (HD 5980) and three in the Large Magellanic Cloud (HD 269006, HD 269333, 
and R 136). The x-axis gives the velocity in km s -1 relative to the Local Standard of Rest; the y-axis shows 
the residual intensity on an expanded scale. 
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Finally, in column 6 of Table 2 we give values of column density 7V(Ti + ), deduced by fitting 
theoretical profiles to the observed absorption lines (Pettini et al. 1990). Details of the profile fits 
(^lsr? velocity dispersion parameter 6, and iV(Ti + )) are collected in Table 3. Since the Ti II lines 
are fully resolved and generally weak, we do not expect significant problems from unrecognized, 
saturated components which, if present, could lead to an underestimate of the total column 
density (Nachman & Hobbs 1973). The errors listed in column 6 of Table 2 reflect only the quoted 
uncertainties in equivalent width. 

Some of our echelle spectra also encompass a resonance line of neutral titanium, Ti I A3341.873, 
with / = 0.1517 (Morton 1991). As expected, this line is not detected since Ti° is a minor 
ionization stage in H I regions; we can place a typical 3cr upper limit of ~ 25% to the contribution 
by Ti° to the total column density of titanium in the H I gas sampled by our observations. 



3. THE SCALE HEIGHT OF INTERSTELLAR TITANIUM 

In order to study the distribution of titanium in the Galactic halo we have combined our 
measurements with those published in the extensive compilations by Edgar & Savage (1989) and 
by Albert et al. (1993). The values of A^(Ti + ) in the compilation by Edgar & Savage have been 
multiplied by a factor of 0.957 to take into account the different /-value of Ti II A3384 adopted by 
these authors. In subsequent figures we have used different symbols to indicate the three sources 
of data, as follows: Edgar &: Savage — circles; Albert et al. — squares; this paper — triangles. 
The combined sample consists of a total of 116 sight-lines; the new observations presented here 
increase the number of objects further than 2 kpc from the Galactic plane from 3 to 12. 

Only one of our 15 stars, HD 220787, is in common with earlier work. Albert et al. (1993) 
report W\(X338A) = 15 ± 5 mA, compared with 26 ± 3 mA measured from our spectrum. 
Accordingly, the value of N(Ti + ) we deduce is ~ 70% higher (Table 2). As this work was being 
completed we were kindly sent a preprint of the recent paper by Albert, Welsh, &: Danly (1994) 
reporting Ti II observations towards six halo stars two of which, HD 195455 and HD 206144, are 
in common with our sample. Again the agreement is not good, with our values of W / a(A3384) 
being ~ 80% higher. These discrepancies are a concern, but we suspect that they may be due to 
the lower signal-to-noise ratio of the spectra by Albert et al. (1993 and 1994). In the present 
analysis we have adopted our own measures for the three stars in question. 

Figure 2 shows the increase in the column density of Ti + , projected along a direction 
perpendicular to the Galactic plane, with distance from the plane | z | for the full sample. It is 
worth noting, in passing, that by considering all the data together, we make the tacit assumptions 
that there are no significant differences in the distribution of titanium above and below the plane, 
nor any trends with Galactic longitude. With a moderately large sample, such as the one under 
consideration here, it is possible to check the validity of these assumptions. Indeed, we find that 
the two sub-samples, at positive and negative Galactic latitudes respectively, are statistically 
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indistinguishable, nor do we see any trends in the abundance of titanium with distance from the 
centre of the Galaxy. Presumably the range of distances probed is too small to show the effects of 
radial abundance gradients (see also Savage & Massa 1985). 



i — i — | 1 — i — | 1 — i — | 1 — i — | 1 — i — 




i i I i i I i i I i i I i i 

1 10 100 1000 10 4 10 5 

|z| (pc) 

Fig. 2. — Column density of Ti + , projected along the z— direction, plotted as a function of distance from 
the plane of the Galaxy. The total sample consists of 116 measurements (or upper limits); different symbols 
denote different compilations of data, as follows: Edgar & Savage (1989) — circles; Albert et al. (1993) — 
squares; this paper — filled triangles. The errors in 7V(Ti+) are those reported by the original authors; the 
uncertainties in the stellar distances are a nominal ±25% . The upper limit at the largest distance refers to 
the sight-line to the Seyfert 2 galaxy NGC 1068; for convenience this point is plotted at | z \ = 80 kpc. The 
increase in 7V(Ti + ) with | z | appears to level off between | z | ~ 1 and 3 kpc. 
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Returning to Figure 2, we see that the increase in iV(Ti + ) with | z | known from earlier studies 
appears to level off between 1 and 3 kpc from the plane, and that the values of iV(Ti + ) measured 
towards the Magellanic Clouds — as well as the upper limit reported by Blades & Morton (1983) 
towards NGC 1068 — are no larger than those for distant stars in the halo. If we assume an 
exponential distribution of the form: 

n 2 (Ti+) = n (Ti+) exp [-| z \/h Ti+ ] (1) 

the column density of Ti + in front of star at height | z | is given by: 

JVsin| 6 | = JV (1 - exp [-| 2? 1/^+]) (2) 

where 

iV = n (Ti+) x h Ti+ (3) 

and b is the Galactic latitude. 

Fitting eq. 2 to the data in Figure 2 and minimizing % 2 in log-log space yielded 

jV (Ti + ) = 1.3 x 10 12 cm" 2 and h Ti + = 1.5 kpc (4) 

For some of the 116 sight-lines considered, only upper limits to iV(Ti + ) are available. We 
treated such upper limits as detections, at a value of iV(Ti + ) a factor of 2 lower than the upper 
limit, and assumed an uncertainty of a factor of 2 (that is, significantly larger than the typical 
errors in iV(Ti+) in cases where the line is detected — see Figure 2). While this procedure is not 
mathematically rigorous, we feel that it is a better approximation than the common approach of 
simply ignoring the non-detections. 

The above fitting procedure also showed that eq. 1 is not an adequate representation of 
the distribution of titanium in the z— direction, since % 2 jv = 14.9, where v is the number of 
degrees of freedom. This was not unexpected, since earlier analyses of this kind have reached 
the same conclusion. It is normally assumed that the large scatter in plots such as Figure 2 is 
due to the intrinsic 'patchiness' of the ISM on the scales sampled, although errors in the stellar 
distances presumably also contribute. We followed Edgar & Savage (1989) in introducing a random 
scattering factor, cr r dm? which was added in quadrature to the individual errors in iV(Ti + ) . We 
then solved for iVo(Ti + ) and h Ti +, forcing y 2 /v to 1 . Table 4 shows the effect on iVo(Ti + ) and 
h Ti + when a r dm is increased from to 0.30 dex; x 2 / u reaches unity when a r dm = 0.29 dex, or 
about a factor of 2. The corresponding fit 

7V (Ti+) = 1.1 x 10 12 cm" 2 and h Ti + = 1.35 kpc (5) 

compared to the values for cr rc jm = 0.00 dex (eq. 4) gives an indication of the uncertainties in 
these quantities. We conclude that, while h Ti + is apparently somewhat greater than the scale 
height of Ca + (h Ca + = 1 kpc; Edgar & Savage 1989), the present data do not support earlier 
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TABLE 4 
Exponential Distribution Models 



All Data (116 sightlines) Galactic Data (111 sightlines) 

a rdm log[7V (Ti+)] h Ti+ X 2 /v log[iVo(Ti+)] h Ti+ X 2 /" 
(dex) (cm -2 ) (pc) (cm -2 ) (pc) 



0.00 


12.12 


1533 


14.94 


12.16 


1770 


15.09 


0.10 


12.06 


1404 


5.12 


12.15 


1753 


5.23 


0.20 


12.04 


1378 


1.86 


12.13 


1702 


1.91 


0.30 


12.03 


1342 


0.94 


12.10 


1604 


0.97 



suggestions that the layer of Ti + in the Galaxy may be several kpc thick. In particular, we note 
that Albert et al. (1994) proposed such an extended distribution of Ti + largely on the basis of 
their observations of one star, HD 123884, of spectral type B8/9Ia/abp apparently located 8.7 kpc 
above the Galactic plane. Albert et al. found Ti II absorption to be exceedingly strong in this 
direction; VF A (A3384) = 172 mA corresponding to A^(Ti+)sin| b \ = 3.9 x 10 12 cm , larger than 
any of the values in our Figure 2. Setting aside the uncertainty in the spectral classification — and 
therefore distance — of the star, the fact that we measure significantly lower values of iV(Ti + ) along 
sightlines to extragalactic sources suggests that either the column density of titanium is unusually 
high towards HD 123884 or the absorption feature identified by Albert et al. is partly of stellar, 
rather than interstellar, origin. 

Returning to our Figure 2, the intrinsic scatter in the distribution of Ti + , cr r dm = 0.3 dex, is 
lower than that of the N(B°) distribution in the same stars (see §4.1 below), cr r d m (H°) = 0.45 dex, 
consistent with the claim by Edgar &: Savage that Ti + is distributed more smoothly than the 
general neutral ISM. 

Finally, we point out that the turn-over in the distribution of Ti + seen in Figure 2 does not 
hinge on the (few) observations of extragalactic sources available at present. If we repeat the 
analysis excluding the four stars in the Magellanic Clouds and the upper limit for NGC 1068, we 
deduce 

AT (Ti+) = 1.4 x 10 12 cm" 2 and h Ti + = 1.8 kpc (6) 

which are only slightly larger than the values appropriate to the full sample (see Table 4). While 
additional observations along sight-lines out of the Galaxy would obviously be very valuable, the 
levelling off of A^(Ti + ) between | z \ = 1 — 3 kpc is already apparent in the spectra of distant halo 
stars. 
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4. THE DEPLETION OF TITANIUM IN THE GALACTIC HALO 

The scale height of Ti + is the largest measured so far for a dominant ion stage in H I gas. 
As discussed in the Introduction, this result is generally interpreted as being due to the reduced 
depletion of refractory elements at the low densities which are presumably found away from the 
disk of the Galaxy. In order to investigate this effect quantitatively, it is necessary to know the 
column density of neutral gas, JV(H°), along the lines of sight studied. 

4.1 Neutral Hydrogen Column Densities 

Values of 7V(H°) are available for 107 out of the 116 directions considered. For nearly all of 
the stars in the compilation by Edgar & Savage (1989) A^(H°) has been deduced from observations 
of the Lyman a absorption line. On the other hand, the values of Af(H°) reported by Albert 
et al. (1993) are from 21 cm emission spectra with an angular resolution of 21 arcmin. These 
measurements may include material beyond the stars, and therefore not associated with the 
observed Ti II, and may also be affected by dumpiness of the ISM within the radio beam. However, 
the corrections required do not appear to be very large, at least for the six stars in the Albert et 
al. survey which have subsequently been measured at Lyman a by Diplas &: Savage 1994b (see 
also Figure 5 of Danly et al. 1992). Specifically, for five of these stars (HD 60848, HD 119608, 
HD 120086, HD 149881 and HD 203664), N(B°) Lya is smaller than 7V(H°) 2 i cm by between « 5 
and « 30%, while for the sixth (HD 214080) N(B ) Lya is « 25% greater than N(B ) 2 i cm (this 
is because Albert et al. attempted to correct for H I beyond the stars by considering only 21 cm 
emission over the velocity range encompassed by the optical absorption lines). In all six cases we 
have adopted A^(H°)i, yct . Lyman a absorption measurements are available for 7 of the 15 stars 
we have observed (see Table 2); for the Galactic sight-lines to 3 of the 4 stars in the Magellanic 
Clouds values of A^(H°)2i CTO have been reported by McGee, Newton, & Morton (1983). Finally 
on this point, molecular hydrogen has been observed in only a subset of the 107 stars; however, 
the molecular fraction of gas in the halo is likely to be small (Savage et al. 1977). 

One way to assess the global properties of our sample of 107 iV(H°) measurements is to 
determine the the vertical distribution of the neutral gas (in a plot of Af(H°)xsin| b | vs. | z | — not 
reproduced here) and compare it with that of the larger data set of 374 sight-lines considered by 
Diplas & Savage (1994a). The best fit to our hydrogen data assuming an exponential distribution 
is obtained for 

7Vo(H°) = 2.3 x 10 20 cm" 2 and h H o = 95 pc. (7) 

This value of No(K°) corresponds to a mean volume density in the disk no(H°) = 0.78 cm -3 ; for 
comparison Diplas &: Savage found no(H°) = 0.366 cm -3 and h H o = 195 pc. Evidently the sample 
of stars for which Ti II measurements are available is less biased against large, dense interstellar 
clouds than the IUE H I survey by Diplas & Savage. 
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4.2 Density -dependent Depletion of Titanium 

Figure 3 shows how the observed ratio of column densities of Ti + and H° varies with distance 
from the plane of the Galaxy. For comparison, the solar abundance is A(Ti)q = 8.5 x 10~ 8 
(Anders & Grevesse 1989), beyond the limit of the plot. The lines drawn are the predictions of 
a simple model of the depletion of Ti discussed below. From Figure 3 it appears that the large 
and highly variable depletion of Ti, well known from earlier surveys, is confined to the disk of 
the Galaxy. The median A(T'i)jsm = 4.2 x 10~ 10 for stars with | z \ < 250 pc corresponds to a 
logarithmic depletion D(Ti)= —2.3, where 

£>(Ti) = log [A(Ti) JSM /A(Ti) Q ] (8) 




I I I I I I I I I I I I I I 

1 10 100 1000 10 4 10 5 

|z| (pc) 

Fig. 3. — Values of Ti abundance measured towards 107 stars at distances | z | from the Galactic plane. 
The symbols have the same meanings as in Figure 2. The four lines drawn give the predictions of a simple 
model for the density-dependence of the gas-phase abundance of Ti, discussed in the text, with power law 
index k = —0.6 (dot-dash line), —0.8 (long dash line), —1.0 (continuous line), and —1.2 (short dash line). 
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Typically, only 0.5% of the titanium is in gaseous form in the disk ISM. In contrast, both 
the average depletion and the scatter are reduced beyond | z | ~ 500 pc; as these sight-lines 
also include material in the disk, the true values in the halo must be even lower. Although the 
number of measurements is not large, there seems to be little residual increase in the integrated 
7V(Ti + )/7V(H°) over sight-lines which extend further than | z | ~ 500 pc; the median value of the 
points at | z \ > 500 pc corresponds to a titanium depletion D(Ti)= —1.5 . 

In Figure 4 the titanium abundance is plotted as a function of the average volume density 
of neutral gas along the line of sight, 7V(H°)/c?, where d is the distance to the star. Similar 
relationships for several other elements have been presented by Jenkins (1987). Again, the lines 
in Figure 4 refer to the model discussed below. The strong trend of increasing A(T\)jsm with 
decreasing < n(H°) > noted by earlier workers can be seen clearly in Figure 4. We deduce a 
gradient of m = —0.70, compared with Jenkins' value of —0.84 . 1 

4.3 Interpretation 

We now attempt to bring together into one coherent picture the information conveyed by the 
data in Figures 2, 3, and 4. We adopt a simple model for the distribution of interstellar titanium 
with distance from the plane of the Galaxy, based on two assumptions: 

(1) The z— profile of H° is represented by the exponential fit 

n 2 (H°) = n (H°) exp[-\ z \/h H o] (9) 

with parameters 

n (H°) = 0.78 cm" 3 ; h H o = 95 pc (10) 
appropriate to our sample (section 4.1 above); and 

(2) the density dependence of the Ti abundance is a power law of the form 

n(Ti+)/n(H°) = A x [n(B°)] k (11) 
where k is negative. Combining eqs. (9) and (11), we obtain: 

n z (Ti+) = A x [n (H°)] fe+1 exp [-| z | x (k + l)/h H o] (12) 
for the gas-phase density of Ti as a function of distance | z | from the plane. 



1 Earlier versions of this plot (e.g. Jenkins 1987; Albert et al. 1993) show an additional point near < n(H°) > = 
0.004 cm -3 and iV(Ti + )/iV(H ) = 4 x 10 -8 , which is clearly separated from the rest of the measurements. We are 
of the opinion that this point, which refers to observations by Stokes (1978) of the sight-line to 77 U Ma, has been 
incorrectly represented. No absorption can be discerned in the original spectrum by Stokes, and the point should 
therefore be considered an upper limit rather than a detection. 



- 18 - 



F~~i — I — in 1 1 1 — I — I I I I I 1 1 1 — I — I I I I I 1 1 1 — I — I I I H 



00 
I 

O 



O 

X 



CD 
I 

o 



o 



o 




_l I 



_l I I I I I I I'll' 



0.01 



0.1 

,0 



10 



<n(H u )> (cm" J ) 



Fig. 4. — Ti abundance plotted as a function of the average density of neutral gas along the 107 sight- 
lines surveyed. The symbols are as in Figure 2; the values of < n(H°) > for the four extragalactic sight- 
lines (plotted with star symbols) assume that the Galactic halo extends to | z \ = 5 kpc. The four lines 
drawn have the same meaning as in Figure 3; they originate at the median value of the Ti abundance 
(A(Ti)isM = 4 x 10~ 10 ) at the base of the halo, where the median gas density is no(H°) = 0.78 cm -3 . 

This model has two free parameters, Aq and the power-law exponent k. By fixing 
Ao x [no(H°)] fc = 4 x 10~ 10 , the mid-plane value of the Ti abundance deduced from the median 
7V(Ti + )/7V(H°) for stars within 100 pc of the plane, and by further imposing the boundary 
condition that A(Ti)ism < A(Ti)q, we can check how different values of k compare with the 
data in Figures 3 and 4 by integrating eq. (12) . The lines drawn in these Figures correspond to 
k = —0.6, —0.8, —1.0, and —1.2 respectively; a reasonable fit is apparently obtained with k ~ — 1 . 
In Figure 5 we show the predictions of the model superposed on the same data as in Figure 2; 
again a value of k ~ —1 is suggested. 

It is important to stress that the comparisons between observations and model predictions 
shown in Figures 3, 4, and 5 are only of a qualitative nature, because of the intrinsic scatter of 
the data and the uncertainties in the choice of model parameters. Nevertheless, the model does 



19 




Fig. 5. — The observed distribution of Ti+ column densities with distance from the Galactic plane is 
compared with those predicted by our simple model of the density-dependence of the Ti abundance. The 
data points are the same as in Figure 2; the four lines drawn are for power-law index k = —0.6 (dot-dash 
line), —0.8 (long dash line), —1.0 (continuous line), and —1.2 (short dash line). 

appear to reproduce broadly the overall trends in Figures 3, 4, and 5, and a rapid increase in Ti 
abundance with decreasing gas density above the plane seems indicated. 

Before discussing the implications of this result, we note that the approach described above 
differs from earlier analyses in one important aspect. Previous workers (e.g. Edgar & Savage 1989) 
have taken the value of k to be the slope in Figure 4. However, that slope measures the ratio of 
column densities integrated along a given distance and is not the same as the index of the power 
law in eq. (11), which is a relationship between the volume densities of Ti + and H° . 



- 20 - 



5. CONCLUSIONS 

We have presented observations of Ti II A3384 absorption towards 15 distant stars in the 
Galactic halo and the Magellanic Clouds. These new data extend earlier surveys of interstellar 
titanium to larger distances from the disk than probed up to now. We report two principal 
findings: 

(1) The scale height of Ti + in the Galaxy is h Ti + = 1.5 ± 0.2 kpc, larger than those of other 
H I species but not extending to several kpc, as had been suspected. The large scale height 
probably reflects the reduced depletion of Ti in low density regions in the halo. 

(2) A simple model in which the density dependence of the gas-phase abundance of Ti is a 
power law with exponent k ~ — 1 reproduces the observations adequately . If this interpretation 
is correct, the volume density of Ti is, on average, independent of the gas density (see eq. 11). 
Apparently, the release of Ti from dust balances overall the falling density of gas away from the 
Galactic plane, until all the titanium is in the gas phase at | z | ~ 1 kpc . This solution is not 
unique. In particular, a faster increase of A(Ti)jsm with decreasing n(H°) is admitted by the data 
if we restrict the maximum gas-phase abundance attainable to be less than the solar value, with 
some of the Ti remaining in solid form even at large distances from the disk. 

It should be relatively straightforward in future to test the validity of our simple model by 
measuring the abundance of Ti directly in individual absorbing clouds which are further than 
| z | ~ 1 kpc. Relevant data are at present available for only a few cases. The clouds at ^lsr = +32 
and +61 km s _1 in R 136 (see Figure 1) are located more than ~ 5 kpc below the plane, if 
their velocities reflect the overall rotation of the Galaxy rather than peculiar motions (Savage & 
de Boer 1981). By comparing our measured A^(Ti + ) = 6.5 x 10 11 and 3.5 x 10 11 cm~ 2 (Table 
3) with N(H°) 21 cm = 2.6 x 10 18 and 7.1 x 10 17 cm~ 2 (McGee, Newton & Morton 1983), we 
deduce D(Ti) = +0.5 and +0.8 respectively (eq. 8). Within the uncertainties introduced by the 
15 arcmin beamwidth of the 21 cm observations by McGee et al., these results are consistent 
with an approximately solar abundance of Ti in these two halo clouds. McGee et al. reached 
a similar conclusion for other refractory elements, such as Al, Si, and Fe, in the halo clouds at 
^lsr = +70 and +120 km s _1 towards HD 36402, also in the Large Magellanic Cloud. As a 
counter-example, in their comprehensive analysis of HST observations of the halo star HD 93521 
Spitzer & Fitzpatrick (1993) found the release of Ti from the grains to be less pronounced than 
that of Fe, with residual depletions .D(Ti) = —0.9 in clouds between «lsr = — 39 and —66 km s _1 . 

Establishing if the solar abundances of refractory elements such as Ti are indeed recovered 
and little dust remains beyond ~ 1 kpc from the plane would provide valuable clues to the 
processes which maintain the interstellar medium in the halo. In particular, the relationship of 
this extended, neutral component to the hot gas producing C IV absorption and emission (see for 
example the review by Savage 1991) is still unclear. While the coronal gas presumably arises in 
a Galactic fountain fed by supernovae (Spitzer 1990; Shapiro & Benjamin 1992; Shull & Slavin 
1994), there may be difficulties in applying the same ideas to explain the distribution of neutral 
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gas within a few kpc of the disk (Houck Sz Bregman 1990). 

It is a pleasure to acknowledge the help by Richard Hunstead and David Mar in obtaining 
some of the data on our behalf, particularly during a brief period in 1990 when one of us (M.P.) 
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very useful in allowing us to complete this project. We thank Geraint Lewis for his assistance 
with the computational aspects of this work, and Barry Welsh for kindly communicating recent 
results in advance of publication. We are very grateful to Blair Savage and Ken Sembach for their 
constructive comments on an earlier version of this paper. 
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